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Abstract 

We present investigations on volume regulation and 
beading shape transitions in PC 12 neurites conducted 
using a flow-chamber technique. By disrupting the cell 
cytoskeleton with specific drugs we investigate the role 
of its individual components in the volume regulation 
response. We find that microtubule disruption increases 
both swelling rate and maximum volume attained, but 
does not affect the ability of the neurite to recover its 
initial volume. In addition, investigation of axonal 
beading -also known as pearling instability- provides 
additional clues on the mechanical state of the neurite. 
We conclude that the initial swelling phase is mechan- 
ically slowed down by microtubules, while the volume 
recovery is driven by passive diffusion of osmolites. Our 
experiments provide a framework to investigate the role 
of cytoskeletal mechanics in volume homeostasis. 

Keywords: Volume regulation, cell mechanics, axon 

beading, cytoskeleton. 

PACS: 87.16.ln, 87.16.dp, 87.16.ad. 



INTRODUCTION 

The ability of living cells to regulate their volume is a 
ubiquitous homeostatic feature in biology Since 
water readily permeates through the cellular membrane, 
alterations in extracellular osmolarity can change the con- 



centration of all cytoskeletal components with severe con- 
sequences for the metabolism. Not surprisingly, one finds 
several mechanisms involved in volume regulation. In 
particular, many eukaryotic cell types display a short-term 
volume regulation response to sudden alterations in exter- 
nal osmolarity, the so-called regulatory volume decrease 
(RVD) and regulatory volume increase (RVI) (i i). It 
is generally accepted that these require the passive dif- 
fusion of ions. In the case of RVD, cell swelling upon 
hypoosmotic shock increases the membrane permeability 
for sodium, which diffuses out of the cell. In turn, the 
decrease of intracellular osmolarity drives water out and 
lowers cell volume In this mechanism water flow 

is driven by a difference in osmotic pressure. It is of- 
ten assumed that hydrostatic pressures are negligible (0), 
with the argument that they would make the membrane 
burst ( 1). However, though the maximal pressures sus- 
tained by lipid bilayers are indeed very low, the situation 
cannot be simply extrapolated to the living cell. The cell 
membrane is connected to the cytoskeleton, the biopoly- 
mer gel spanning across the cell interior (5-7). Since the 
cytoskeleton is viscoelastic and contractile (8-1 1) (also in 
PC12 neurites (Il2h). it may provide a mechanical memory 
of the initial state as well as a driving force for volume 
relaxation (il3i-t20i). Though it has been argued that the 
cytoskeleton is too weak (typical moduli are up to 10 kPa 
(fi 10)) to sustain osmotic pressures (up to '--'MPa) Q), 
neurons under hypoosmotic shock sustain strong pres- 
sures for several hours (21) and swelling of erythrocytes 
increases and approach perfect osmometer behaviour af- 
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ter disruption of the spectrin-actin cytoskeleton (22). This 
is indeed compelling evidence for a mechanical role, but 
in many systems the cytoskeleton seems to additionally 
play a signalling one. In particular, biochemical interac- 
tions between actin filaments and ion channels may cou- 
ple strain of the actin cortex to changes in channel activity 
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231) . The fact that the actin cortex is disrupted when 



hypoosmotic swelling begins (|24|. |25|) seems to be due 
to an influx of Ca^+ (26) through mechanosensitive ion 
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channels activated by membrane stretching (0, ill)- To Figure 1: A schematic diagram of the flow-chamber set-up. 



clarify the role of the cytoskeleton one must discern be- 
tween pure mechanical and mechanosensing responses, a 
difficult task requiring direct measurements of membrane 
tension. 

In this work, we study volume regulation in neurites, 
axon-like cylindrical protrusions extended by PC 12 cells 
(l29b . structurally very similar to the axons produced by 
neurons in culture (l6i[30h. Neurites furnish a hitherto un- 
explored, yet attractive model system to investigate the 
role of mechanical tension in volume regulation. Their 
simple cylindrical geometry and low amount of invagi- 
nated membrane allows better volume and area calcula- 
tions from images. They also have a well-defined, highly 
organized cytoskeletal structure similar to that of axons: 
a central array of longitudinally arranged microtubules 
interconnected by microtubule binding proteins and sur- 
rounded by an actin cortex. Moreover, an exceptional 
feature of axons and neurites is their ability to undergo 
sudden shape transformations in response to an applied 
mechanical tension (2^ 31, 32), an instability known as 
axonal beading to biologists and pearling instability to 
physicists. We show that axonal swelling increases mem- 
brane tension and that microtubules slow down the vol- 
ume response. We conclude that frictional forces in the 
cytoskeleton play an important role in axonal volume reg- 
ulation. 

EXPERIMENTAL SETUP AND METH- 
ODS 

Flow chamber 

The experiments have been carried out using a flow- 
chamber technique. A schematic of the set-up is shown in 
Fig. [T] A stainless-steel block and two coverslips are used 
to form a 10 X 5 X 1 mm'^ chamber. Cells adhere on the bot- 
tom cover-slip. One duct of the chamber is connected to 
a peristaltic pump by means of long, soft silicone-rubber 
tubing that minimizes pressure fluctuations arising from 
the pump. The other duct is connected to a 3-way-valve 



to select between two different media that can be pumped 
into the chamber. Stainless-steel tubes that are 1 mm in 
diameter connect the 3-way-valve to the two reservoirs 
where the media are stored. The chamber is intentionally 
made small to ensure a quick switching from one medium 
to the other at low flow rates, in the range of 2— 4/il/s. The 
switching of concentrations inside the chamber was stud- 
ied by adding an absorbing dye to one of the solutions and 
monitoring the variation in the transmitted intensity with 
time. The concentration reaches 90% of its final value 
in about 10 s. The chamber, the 3-way-valve, and the 
stainless-steel tubes are placed inside an aluminium block 
with good thermal contact between each other and a wa- 
ter bath maintains the temperature of the block at the de- 
sired value. The continuous flow of pre-warmed medium 
keeps the chamber at constant temperature despite some 
heat loss through the cover-slips. 

Cell culture 

PC 12 cells are from the German Collection of Microor- 



ganisms and Cell Lines (DSMZ) (l33h . They are plated 
on collagen coated slides and cultured in RPMI-1640 
medium (Gibco) with 10% fetal bovine serum and 5% 
horse serum in presence of nerve growth factor (NGF) 
(Sigma- Aldrich Chemie, Munich, Germany) for 4-5 days 
(jSJ). Such young neurites are known not to have in- 
termediate filaments (130!). Slides are coated with 3- 
aminopropyl triethoxysilane (Sigma). The silanised slides 
are covered with about 100 /il of 10% rat-tail collagen 
(Sigma- Aldrich) dissolved in a 70% ethanol - 30% water 
solution and let dry overnight. 

Experimental procedure 

Prior to an experiment, the slide with the adherent cells 
is transferred to the flow-chamber. Cells are allowed to 
stabilize for about 5 min by circulating the experiment 
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medium (normal medium with addition of 25 mM HEPES 
buffer (Invitrogen, Darmstadt, Germany)). Experiments 
are performed by switching from the normal experiment 
medium to experiment medium diluted with deionized 
water. The response of the cells to the sudden switch in 
the external concentration is observed with an Axiovert 
135 microscope (Zeiss, Oberkochen, Germany) config- 
ured for phase-contrast imaging. 

Image analysis 

The volume and area of the neurites are analyzed from the 
recordings using a home-made edge detection program. 
Edge detection using a threshold for intensity is unreli- 
able due to the "halo effect" present in phase-contrast im- 
ages and also due to the dependence on the illumination 
intensity. To avoid such complications, edge pixels are 
recognized along the neurite by detecting the local max- 
ima in the gradient of intensity across the neurite. Af- 
ter edge detection the neurite volume and surface area 
are computed assuming axial symmetry for the neurite 
shape. Axial symmetry is verified by comparing the two 
detected boundaries of several neurites and is found to be 
a good approximation for straight neurites which are at- 
tached only at the two extremities. Only such neurites 
were selected for the experiments. 

Drug-induced cytoskeletal perturbation 

Experiments were performed in presence of cytoskeleton 
disrupting drugs in order to study the role of its individ- 
ual components. A complication in these experiments 
arises due to the neurites becoming fragile or losing their 
cylindrical geometry on cytoskeletal perturbation. This 
precluded experiments with the actin disrupting drug La- 
trunculin, which induces detachment of the growth cone. 
In the case of the microtubule disrupting drug Nocoda- 
zole (NOC) ("ii"), which induces shape irregularities after 
approximately 10 min exposure, we let the drug act for 
5 min before performing the hypoosmotic shock. NOC 
concentration was 10 /Ug/ml throughout. In contrast, the 
myosin-blocking drug Blebbistatin (36) did not signifi- 
cantly alter neurite shape. Thus, to ensure its effect we in- 
cubated neurites at 37° for 1 hour at a high concentration 
(50 jjM) before transferring them to the experiment cham- 
ber and performing the experiment at a lower concentra- 
tion of 20 fjM. For all experiments, drugs were present 
both in the normal and in the diluted medium. Since all 
drugs are dissolved in dimethylsulfoxide (DMSO), which 
itself has effects on water and ion channels (37), we per- 
formed control experiments in presence of DMSO 0.5%, 



equal to the highest DMSO concentration in any of the 
drug experiments. 

EXPERIMENTAL RESULTS 
Volume dynamics 

We begin all experiments with a hypoosmotic shock 
imposed by switching the cell culture medium flowing 
through the chamber from normal medium with a total 
solute concentration Co ^ 300 mM to a lower external 
value Ce (the continuous flow ensures constant external 
concentration at all times). In the following, for an intu- 
itive measure of the shock strength we will normalize the 
external osmolarity by its initial value and denote it by 
C = Cs/Cq. Figure [2] shows typical responses for three 
different values of C at 36°C. For weak shocks (C = 0.8) 
the neurite volume increases from its initial volume Vq un- 
til it reaches a maximum steady value Vmax- No recovery 
is observed for tens of minutes. For intermediate shocks, 
C = 0.7, the volume increases at a roughly constant rate 
Vq initially until it reaches a maximum value Vmax- Sub- 
sequently, the volume recovers almost back to Vq with a 
typical regulatory volume decrease response (RVD). The 
volume recovery is roughly exponential with a character- 
istic time T. The recovery time r is strongly temperature- 
dependent (shown later), but does not follow a clear trend 
with neurite diameter. For strong shocks, C = 0.5, the 
recovery is faster and there is often a remarkable "un- 
dershoot", where the volume goes below its initial value. 
Moreover, during the swelling phase the neurite under- 
goes a sudden and transient shape transformation from its 
normal cylindrical geometry to a periodically modulated 
one. As discussed later, this Rayleigh-like "pearling" in- 
stability or beading indicates an increase in membrane 
tension upon neurite swelling. This shape transformation, 
which will be discussed later, does not affect the volume 
response curve in any measurable way. This fact allow us 
to discuss the volume dynamics and shape dynamics in 
that order. 

Once the volume stabilizes to the lower external osmo- 
larity Ce (within about ten minutes), we perform a hy- 
perosmotic shock by switching back the original medium 
with concentration Cq. The neurite shrinks, reaches a 
minimum volume Vmin and then comes back to its ini- 
tial volume in an RVI response as shown in Fig. |2] No 
beading shape transformation is observed in this case. 

In the following, we separately address the volume reg- 
ulation response and the peristaltic modulation. We study 
them as a function of osmotic shock strength, tempera- 
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ture, and in the presence of cytoskeleton-disrupting drugs. 

Nonlinear swelling response. Figure [3] A shows the 
initial rate of change of volume divided by the initial area, 
Vq/Aq, as a function of the initial osmotic pressure dif- 
ference Alio = RTCo{l — C), for temperatures 33- 
36°C. Assuming zero hydrostatic pressure and neglecting 
changes in internal osmolarity gives the following expres- 
sion for the initial swelling rate: 

yo = ^o^pAno, (1) 

where Lp is the hydraulic permeability. Surprisingly, 
from the average Vq value at C = 0.7 we obtain an os- 
motic permeability Pj = RT6w Lp ~ 1.4/im/s (where 
6w is the molar density of water), which is about two 
orders of magnitude lower than the literature values for 
lipidic membranes and most biological cells ^ [381-0) 
even after blockage of water channels (37). Moreover, as 
can be seen from Fig.|3]A, the expected linear dependence 
is contradicted by the strong nonlinear response observed 
for both hypo- and hyperosmotic shocks. 

We now turn to the maximum (minimum) volume at- 
tained in a hypoosmotic (hyperosmotic) shock. The max- 
imum (minimum) volume Vmax{Vmin) is to a good ap- 
proximation proportional to the initial volume Vq and 
does not depend significantly on the temperature (data 
not shown (41)). Thus we look at the relative maximum 
volume Vmax/yo in order to minimize the effect of neu- 
rite diameter. As shown in Fig. [3]B, the maximum vol- 
ume increases nonlinearly with the initial osmotic pres- 
sure difference Alio. The data is contrasted to the per- 
fect osmometer equation, corresponding to a constant to- 
tal amount of internal osmolites and zero hydrostatic pres- 
sure: 

t7 tit ~ ' 

^max •'dead 

where the dead volume Vdead represents non-aqueous in- 
ternal volume. Mammalian cells have on the average a 
cytosolic protein concentration of ~ 20% (5J. Accord- 
ing to electron microscopy studies the non-cytosolic 
volume of PC 12 neurite is comprised mostly of micro- 
tubules and organelles and amounts to Vdead — 25% Vq. 
The black region in Fig. [3] B corresponds to y > Vq/C 
for hypoosmotic shocks and to y < Vq/C for hyperos- 
motic shocks; penetrating this region would require work 
against the osmotic gradient. As expected for a passive 
response, the data lies outside them. For mild shocks, 
C > 0.7 and (hyperosmotic) C < 1.4, neurites swell 
as much as perfect osmometers a dead volume of 25%. 
Therefore, any ion leakage or hydrostatic pressure are 
negligible during the swelling phase. However, at strong 
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Figure 2: A: Swelling and recovery after a hypoosmotic shock 
performed at time t — Q. Normalised volume V/Vq as a func- 
tion of titne for different dilutions: C = 0.5, (circles), C ~ 0.7 
(solid line), C — 0.8 (dashed line). Temperature is in all cases 
36° C. Each curve corresponds to a different neurite. The max- 
imum volume attained is Vmax, ond the minimum is Vend- No- 
tice the strong undershoot of the volume at C = 0.5, and the 
absence of recovery atC — 0.8. These are all general trends of 
the volume response. B: The full observed response when the 
neurite is subjected to a hypoosmotic shock and a subsequent 
hyperosmotic shock. At first the hypoosmotic shock Cq Cg 
is applied and the neurite is allowed to swell and relax (RVD). 
After the relaxation phase the hyperosmotic shock Ce — ^ Co 
is applied and the neurite shrinks and recovers the initial vol- 
ume (RVI). Such Co ^ Ce ^ Co cycles can be repeated up 
to five times in a given neurite before the ends detach (data not 
shown). 
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Figure 3: A: Initial swelling speed Vq/Aq as a function of 
the initial osmotic pressure difference Alio for temperatures 
33 — 36°C. The shaded region is a guide to the eye, a spline 
going through the averages within a standard deviation. No- 
tice the nonlinear response to strong hypo- as well as hyperos- 
motic shocks. B: Maximum change in relative volume Vai/Vq 
as a function of the inverse external concentration 1/C, for 
all temperatures (here Vm stands for Vmax or Vmin os the 
case may be). The grey region is a guide to the eye. The 
dashed line corresponds to a perfect osmometer with dead vol- 
ume Vdead — 0.25Vo, and the black regions correspond to neg- 
ative dead volumes. 



shocks (hypoosmotic C = 0.5 or hyperosmotic C = 2), 
neurites swell significantly less than a perfect osmome- 
ter would. Thus, upon strong osmotic perturbations either 
ions leak or a sustained hydrostatic pressure develops. 

Volume regulation under cytoskeleton disruption. 

The axonal cytoskeleton may be expected to contribute 
to the volume response in several ways. As discussed in 
the introduction, a mechanical as well as a signalling role 
is conceivable. In order to assess the role of individual 
components of the cytoskeleton, we treat neurites with 
the myosin blocking drug Blebbistatin (BLE) (1361) and 
the microtubule disrupting drug Nocodazole (NOC) (^ 35|) . 
Since these are all diluted in dimethylsulfoxide (DMSO), 
a compound known to alter ion channels, we also perform 
control experiments in presence of DMSO. 

Fig- HI A shows typical responses. In presence of Noco- 
dazole the initial swelling rate for strong shocks C = 0.5, 
increases markedly, but it barely changes for C = 0.7 
(see Fig.|4]B). With Blebbistatin we observe a weaker but 
still significant effect. For both drugs, the relationship be- 
tween swelling rate and initial osmotic pressure difference 
approaches the naively expected linear dependence given 
by Eq.[I] 

Nocodazole induced disruption of microtubules also 
has a strong effect on the maximum volume Vmax attained 
after a strong shock of C = 0.5, as shown in Fig. |4]C. 
Neither BLE-treatment nor DMSO alone have a signifi- 
cant effect on the maximum volume. For mild dilutions, 
C = 0.7, nocodazole has no effect, consistent with the 
fact that neurites swell like a perfect osmometer with 25% 
dead volume. 

Importantly, the cytoskeleton disrupting drugs do not 
affect the ability of the neurite to perform RVD. As can 
be seen in Fig. |4] A, the volume fully relaxes back to its 
initial value. Further evidence can be found in the Sup- 
plementary Material. 

The effect of temperature: Arrhenius behaviour. 

We now address the influence of temperature in the dy- 
namics of volume regulation. For simplicity, we describe 
the recovery phase by fitting single exponentials. As 
shown in Fig. |5l lowering the temperature from 35 to 
15°C slows down the volume dynamics by an order of 
magnitude. The dashed line is an Arrhenius-like equation 

AG 

T (X 1/k (X CRT , where k depends on temperature as the 
rate constant of a thermally-activated process. This yields 
an activation energy AG ~ 30 kT, a typical order of mag- 
nitude for biological processes (|42j). The well-defined Ar- 
rhenius trend is consistent with the idea that ion channels 
are responsible for RVD. 
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Figure 4: A: Comparison of different drug treatments. Each 
curve is a different neurite. All experiments were performed at 
temperature T — 33°C and dilution C — 0.5. Grey shaded 
line: DMSO control. Dashed line: BLE. Black solid line: 
NOC. B: Swelling speed Vq/Aq as a function of the initial os- 
motic pressure difference AHq, for temperatures 33 — 36°C. D.- 
DMSO control (n = 5 for C = 0.5, n = Ufor C = 0.7). B: 
blebbistatin (n = 7, n = 13j. N.- nocodazol (n = 8, n = 14j. 
The grey region is a guide to the eye, corresponding to the ex- 
periments without drugs shown in Fig. \3\A. Data for all drugs is 
shown as arithmetic mean ±S. E. C: Maximum relative volume 
Vm /Vq as a function of the osmotic pressure difference AH 
in presence of cytoskeleton-disrupting drugs. The grey region 
corresponds to the data without drugs in Fig. \3\B. 
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Figure 5: Recovery time as a function of inverse temperature 
1/T. The dashed line is a least-squares fit to an Arrhenius form 



Beading 

When neurites are subjected to strong hypoosmotic 
shocks they undergo a shape transformation by develop- 
ing a periodic array of swellings akin to beading of axons 
(1431) . We observe this in chick dorsal root ganglia (DRG) 
neurons as well as PC 12 neurites Fig. [6l This peristaltic 
deformation, resembles that formed in nerves subjected 
to induced stretch injuries and that observed 

after traumatic injuries to the brain (43). The dynamics 
of bead formation and the mechanism has been investi- 



gated recently using the osmotic shock technique (|2Clh and 



shown to be similar to the pearling instability observed in 
synthetic membrane tubes under tension ^A^. Here, 
we provide direct evidence correlating neurite tension and 
bead formation and describe the role of cytoskeletal com- 
ponents in this process. 

We begin by listing the main features of osmotic shock 
induced beading. For a given radius of the neurite, there 
is a critical hypoosmotic shock below which the shape 
remains cylindrical during the entire volume evolution, 
and above which a transient peristaltic modulation is ob- 
served. This is about C = 0.5 at 37°C and about C = 0.7 
at 25°C for an initial neurite radius of 0.7 ^m. Beading 
and recovery cycles (for mild shocks) can be repeated up 
to five times in the same neurite, after which neurites tend 
to detach from the substrate. Also transport of organelles 
can be observed during and after beading and recovery. 
This suggests that mild osmotic shocks leading to bead- 
ing and recovery causes no permanent damage to the neu- 
rites. On the contrary, strong shocks {C > 0.5) often lead 
to strongly modulated beaded shapes with no recovery for 
several tens of minutes. Extreme shocks cause neurites to 
burst with organelles spewing out, corroborating a build 
up of membrane tension. Finally, no beading is observed 
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Figure 6: Top: Swelling-induced pearling instability, a: 
Chick-embryo neurons, b: PC12 neurites at a higher magni- 
fication. Bottom: Pearling instability. Relative volume V/Vq, 
area A/Aq and shape parameter ^JV/Vq / {A / Aq) as a func- 
tion of time. Imposing a fast hypoosmotic shock at dilution 
C = 0.5 and temperature 33°C increases membrane tension 
and leads to a shape instability. The instability vanishes before 
the volume recovers. 




Figure 7: Neurite deformations induced by the combined ef- 
fect of a constant laminar flow perpendicular to the neurite 
and different osmotic conditions. The temperature is 25° C and 
Ce — O.TCq. The sequence of events is as follows: laminar 
flow is started at t = 20 s and the neurite allowed to undergo 
viscoelastic relaxation and attain a steady state, a hypoosmotic 
shock is applied at t = 320 s by introducing dilute medium, 
the neurite is allowed to undergo volume and shape relaxation, 
normal medium is reintroduced after the volume has relaxed 
(hyperosmotic shock) at t = 1200 s, the neurite is again al- 
lowed to undergo volume relaxation. The flow rate is constant 
throughout the experiment. The dashed lines are for compari- 
son of curvature. The larger the average curvature of the cate- 
nary, the smaller the tension in the neurite ^4% . 



during hyperosmotic shocks, irrespective of the magni- 
tude. 

Axonal tension causes beading. The flow chamber 
technique has been used recently to quantitatively study 
the evolution of neurite tension in response to drag forces 
and to demonstrate active contractile responses in neurites 
(|47h . This is performed by imposing a constant, laminar 
flow perpendicular to the neurite generating a drag force. 
It has been shown that PC 12 neurites as well as axons 



are viscoelastic (Il2h and respond to a stretching force via 
a relaxation process with a characteristic timescale. Un- 
der the influence of a flow induced drag force the neu- 
rite takes the form of a catenary and elongates until it 
reaches a final equilibrium shape (0)(seeFig.|7]t=3OOs). 
The resulting strain gives a measure of the tension in the 
axon. Fig. |7] shows snapshots of the neurite taken at var- 
ious stages during an osmotic shock experiment. Under 
the influence of the flow the neurite attains an equilibrium 
catenary shape within about a minute. When the flow is 
switched from the normal medium to a diluted one, keep- 
ing the flow rate the same, the neurite shortens until it 



becomes almost straight. The beading sets in during this 
straightening phase and the maximum beading amplitude 
occurs when the neurite is straight (t=380s in Fig.©. Sub- 
sequently there is a relaxation process in which the bead- 
ing amplitude begins to decay and the flow-induced cur- 
vature increases. When a hyperosmotic shock is apphed 
after the first volume relaxation process (t= 1200s in Fig. 
IT) the neurite curvature begins to increase as the neurite 
shrinks and subsequently decreases to almost its normal 
value with a timescale very much comparable to the vol- 
ume recovery time. These observations clearly demon- 
strates that the neurite tension and the volume changes 
induced by osmotic treatments are correlated. 

Beading mechanism. The physical mechanism for 
bead formation is as follows (l20r) . After applying an os- 
motic shock the neurite volume increases as discussed 
above. This results in a corresponding expansion of the 
measured area A{t) as shown in Fig. [6l This stretch- 
ing of the membrane causes an increase in the tension 
of the outer membrane and hence an increase in the sur- 
face energy, Fg = CT{A{t) — Aof' / Aq. The increase in 
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volume also costs bulk elastic energy due to the deforma- 
tions caused in the cytoskeletal network, which we will 
consider below later. 

It can be shown (46) that at any instant t, for a volume 
V{t), a peristaltically modulated shape with 

r{t, z) = rQ + e{t) sin{qz) 

has a lower average surface area compared to a cylinder 
with the same instantaneous volume V{t), provided the 
wavelength A > 2iTrQ. For small amplitudes, the relative 
area gain can be obtained as 

5S/S = e{f - l)/(4r2) , 

where S is the surface area and q = lirr/X. This can 
be seen in Fig. [6l The ratio ■\/V{t)/A{t) (note: the av- 
erage volume and area are computed for unit length) in- 
creases when the peristaltic mode grows. This quantity is 
a constant for all cylinders irrespective of radius. An in- 
crease in 'syV{t)/A{t) from this constant value indicates 
a decrease in area for the peristaltic shape as compared 
to a cylinder with identical volume. In other words, the 
neurite is able to reduce its interfacial energy by adopt- 
ing a peristaltic shape instead of a cylindrical one. As 
mentioned earlier any deformation costs bulk elastic en- 
ergy due to the cytoskeletal elasticity. Unlike the surface 
free energy, the bulk elastic free energy is always positive 
for peristaltic modes compared to a cylinder of same vol- 
ume (it costs bulk energy to expand as well as compress 
regions along the neurite). Therefore, there is a compe- 
tition between the surface energy and bulk energy which 
determines the preferred shape, giving rise to a critical 
tension. A rough expression for the critical tension can be 
obtained by comparing the average bulk energy per unit 
length of a cylinder Er'^, where E is the elastic modulus 
of the cytoskeleton, and the corresponding surface energy 
avc- This gives a critical tension Uc = Er^, above which 
the surface contribution dominates. As seen in experi- 
ments, the critical tension needed to destabilise the cylin- 
der increases with neurite radius. Analogous shape trans- 
formations observed in membrane tubes, triggered by ap- 
plication of laser tweezers (,45^ 46.). a nd in cell protrusions 
after treatment with latrunculins (|48h have been analysed 
in a similar fashion. 

The above mentioned expression for relative gain in 
area shows that reduction in area increases with wave- 
length of the perturbation, giving a maximum area gain 
for A = L/2, L the length of the neurite. Clearly, the 
observed value of A is much shorter than and independent 
of the neurite length, but increases linearly with radius. A 



simple minded argument for this is as follows. Any per- 
turbation with a wavelength A generates as pressure dif- 
ference given by the Laplace law, the pressure difference 
p balanced by the membrane tension a. In the case of the 
axisymmetric modes there are two principal curvatures: 
that with a curvature radius equal to the neurite radius 
r{z), and that corresponding to the peristaltic deformation 
along the z-axis, with a curvature radius l/d1r{z). For 
large wavelengths, A » e, the latter can be ignored. The 
pressure difference in between the crest and the trough of 
a peristaltic mode can be written in the small amplitude 
limit (e <^ r) as 

= a{l/rcrest - l/rtrough) ^ -cre/r'^ . 

As indicated by the negative sign, this is an unstable flow: 
the pressure difference will drive water from the troughs 
into the crests, thereby increasing the perturbation ampli- 
tude e. The flow rate and hence the growth of a given 
mode is in general proportional to the driving force 

Ap/X ~ ae/ Xr'^ . 

Shorter wavelength modes are thus faster to grow than 
longer wavelength ones. However, for very short wave- 
lengths the curvature along the neurite axis becomes im- 
portant. The respective pressure difference is given by 



Ap 



ad;r[z) 



ae/X . In contrast to the previous 



case, here the cylindrical state is stable: water flows from 
the crests into the troughs (which have a "negative" pres- 
sure) making the perturbation vanish. Moreover, the de- 
pendence of the flow rate on the wavelength is more pro- 
nounced, 

Ap/X ~ ae/X^ . 

Therefore for very small wavelengths the stabilizing flow 
driven by the curvature along z always dominates, and 
there is no instability anymore. We conclude that there is 
a fastest mode at an intermediate, non-zero wavelength. 
This can be shown to be A ~ 9.2 r (l49l) . In the very dawn 
of linear stability analysis, this argument was applied to 
the case of water jets by Lord Rayleigh (I50l. l5lh. 



THEORY 

We model volume regulation in neurites taking into ac- 
count both mechanical and osmotic driving forces. The 
neurite is characterised by its volume, V{t), and the in- 
ternal amounts of ionic species, ni{t). For simplicity we 
consider only potassium and chloride, by large the most 
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important osmolites inside the cell. Since due to elec- 
troneutrality the flows of Cl~ and K+ are coupled, 

"-C1 = riK , 

one of the two concentrations can be eliminated. The two 
variables of our model are the adimensional quantities 

V = V{t)/Vo (3) 
Af = nci(t)/nci(0) . (4) 

The flow of water through the membrane 111) is gi 



nven 



by 



V 



ALp{AU - Ap) , (5) 

where A is the neurite area, Lp the hydraulic permeabil- 
ity, Ap the hydrostatic pressure, and the osmotic pressure 
difference is given by 



An~i?r^(ni/y-c| 



in terms of external c^^ and internal osmoUte concentra- 
tions Ui/V. The typical relaxation time for the volume 
is 

^0 



TV 



ALpUo 



where Ho = RTCq ~ 700 kPa the osmotic pressure of 
normal medium. Taking 10~^^ m Pa^^ s^^ for the per- 
meability (from Fig.|3]A) we obtain Ty = 50 s. 

We model ion movement with a passive K+-C1^ co- 
transport following Ref. (3), neglecting the effect of ion 
pumps since these are not relevant for the short-term RVD 
response. Chloride flux is given by the difference in 
chemical potentials. 



raci = -ARTG\og 



/ ncinK 



1 



(6) 



We model the RVD response following the standard as- 
sumption of a volume-dependent permeability ^. The 
permeability G must be zero at the initial volume Vq and 
non-zero at a significant departure. We decompose 

G{V) = \G\ g{V) , 

into a typical order of magnitude |G| and an adimensional 
function g{V). For the chloride relaxation time 



Vr 



TCl 



^Cl 



ART\G\ 



we expect a value of 10-100 s taking typical values for 
ion conductivities (1). The volume-dependent permeabil- 
ity function g{V) jumps from zero to a finite value at a 



critical threshold volume; otherwise its form is not known 
and different expressions have been used in the literature 
(0, [l7l) . In our experiments the critical volume is about 
1.2 Vq, as can be seen in Fig. |2] We obtain good results 
with the simple expression 

g = e{V/1.2)V, 

where 9{) is the step function. 

Finally, assuming external ion concentrations to jump 
instantaneously at time t = from c?^(0) to Cc?^(0) the 
full system can be written as 



V 



J_ / 1 + 20(AA- 1) Ap 

^ V V no 



9{V) 



TCI 

with the par ameters 

</> = 
a = 

/? = 



log 



a 



AA2 + (l 



C 



/3C2V2 



An I Jn 



•-CI / "^K 
ex ex / in in 



(7) 

(8) 

(9) 

(10) 

(11) 



Taking the following physiological values for the ion con- 
centrations (Jj) 





= 80mM 


Jn 


= 150 mM 


Jn 


= 15 mM 


cci 


= 80mM 


^ex 


= 6mM 


^ex 
CNa 


= 110 mM 



the adimensional constants become (\> = 0.25, /3 ~ 0.05, 
and a ~ 0.5. 

Zero pressure model. The standard assumption in 
modelling volume regulation is to set the hydrostatic pres- 
sure to zero (l3iD, 

Ap = . 

With this ansatz we can solve the equations for the volume 
response V(t). As shown in Fig.Hl the model reproduces 
the timescales for swelling and recovery as well as the 
maximum volume. Interestingly it reproduces the "un- 
dershoot" at C = 0.5 (see Fig. |2]l. However, the depen- 
dence of the swelling rate on osmotic pressure is clearly 
off. The model gives a linear response -a feature intrinsic 
to any volume-dependent permeability- whereas experi- 
ments show a nonlinear one (Fig. [3] A). 
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V/Vo 




time (s) 



time (s) 



Figure 8: A: zero pressure model for dilutions 80% (dotted 
line), 70% (dashed line) and 50% (solid line). B: adding non- 
linear friction as given by Eq.\T2\reproduces the trend during 
the swelling phase. Compare to Figs. \ 



Viscoelastic model. As shown by our experiments 
combining hypoosmotic shocks and drag forces (Fig. O, 
neurites develop tension while swelling. This suggests 
that the nonlinear response observed in experiments may 
come from a nonlinear mechanical response, a ubiquitous 
feature of cells with structured cytoskeleton (.52<) . We 
therefore assume that the pressure difference Ap is not 
zero, but depends on the swelling rate as expected for a 
viscoelastic element. To reproduce Fig. |3]A, we assume 
that neurites swell at zero pressure below a critical rate 
but encounter internal friction for rates larger than a criti- 
cal value Jl. This can be written as 







r/sign(V) {\V\/^-l) 



for 
for 



1^1 



(12) 



where is the maximal rate of change of volume at which 
the neurite can swell without tension, corresponding to 
about 5 10~^ in our experiments, and rj is the friction 
scale, about 2 MPa. 

As shown in Fig. [H with this ansatz the model re- 
produces the essentially constant swelling rate observed 
in experiments. The nonlinear friction slows down both 
swelling and recovery. Interestingly, the peculiar triangu- 
lar shape of of the curves at C = 0.5 is much closer to 



that of experiments (compare to Fig.|2]). 

DISCUSSION 

Microtubules mechanically slow down swelling. The 

initial response of neurites to hypo- as well as hyperos- 
motic shocks is a strongly nonlinear function of the ex- 
ternal osmotic pressure. Remarkably, this response be- 
comes much simpler after microtubule disruption: neu- 
rites behave as perfect-osmometers, and the swelling rate 
increases linearly with osmotic shock strength. The sim- 
plest explanation is that microtubules slow down the ini- 
tial volume change by mechanically opposing the hydro- 
static pressure difference. This is not exactly an elastic 
response but rather a (nonlinear) viscous one which de- 
pends crucially on the swelling rate. This mechanical role 
of microtubules may come as a surprise, as they are in 
general irrelevant for the (passive) mechanical response 
of cells. In neurites, however, their structure is differ- 
ent: they are arranged in bundles interconnected by mi- 
crotubule binding proteins (MBPs). Our results indicate 
that microtubules are firmly connected to the membrane, 
either directly by MBPs or through the actin cortex. In- 
terestingly, studies of volume regulation on round PC 12 
cells without neurites give a different picture. Whereas 
disruption of the actin cytoskeleton by Cytochalasin B has 
been reported to increase KCl efflux and diminish extent 
of swelling / isi Isih . disruption of microtubules has no 
effect on RVD (l25b . This is consistent with our conclu- 
sion of a mechanical role for micro tubuli, since these are 
organised very differently in round cells and neurites: in 
the latter their bundle structure provides a rigid scaffold 
which can oppose swelling. 

The mechanical response of the neurite, as described 
by Eq. [T2l se ems to be similar to that of adhering fibrob- 
lasts (1521 p4h . At slow strain rates, forces are low; above 
a critical strain rate, friction increases and the force be- 
comes much stronger. From our results, writing the crit- 
ical rate as a strain we get Vb/(''o^o) ~ SlO^'^s^^. In 
single fibroblasts, this is indeed the order of magnitude 
of the critical strain rate where frictional forces increase 



(|52|) . which in turn agrees with the timescale of active pro- 
cesses (0, [ill) . This suggests that changes in cell shape 
-including swelling- take place at the rate allowed by 
spontaneous unbinding of cytoskeletal crosslinks, while 
faster changes are slowed down by friction between cy- 
toskeletal elements. 

RVD is osmotically driven. Our results indicate that 
mechanical tension does not provide the driving force 
for volume recovery. The pearling modulation vanishes 
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well before the volume recovers, indicating zero mem- 
brane tension during RVD. Moreover, in presence of 
cytoskeleton-perturbing drugs the volume recovery time r 
does not increase, but becomes slightly shorter. This indi- 
cates that volume recovery takes place via extrusion of os- 
molites Indeed, our simple model based on two ionic 
species suffices to explain the essential features of the vol- 
ume evolution curves. It captures the broad response and 
even the tendency to "undershoot" at large dilutions. 

RVD follows an Arrhenius trend as a function of 
temperature. The timescales involved in the volume re- 
sponses are temperature dependent. The relaxation time 
is about an order of magnitude more temperature depen- 
dent than the swelling/shrinking rate, showing an expo- 
nential reduction as the temperature is increased. Thus the 
ability of the cell to respond to and negate perturbations 
to its volume improves drastically as the temperature ap- 
proaches the normal physiological values. Due to this rea- 
son the neurites are able to withstand strong hypoosmotic- 
shocks at higher temperatures when the same shock 
would have made them rupture at a slightly lower tem- 
perature. 

A dynamic picture of stretch injury. Similar shape 
transformations have been observed in nerves under the 



name of "beading" as a response to stretch injury (|3ll. l32h. 
Interestingly, electron microscopy observations of the ul- 
trastructure of stretch-beaded nerves show that micro- 
tubules are splayed out in the beads ^4^, consistent with 
our picture of a mechanical connection between micro- 
tubules and the membrane. The shapes of beaded nerves 
have been interpreted as equilibrium shapes with a con- 
stant curvature (31); however, our results offer an alter- 
native explanation, namely that the shape is given by the 
fastest growing mode at the time of the increase in ten- 
sion. If so, the shape of beaded nerves would not simply 
follow from structural properties (as expected for an equi- 
librium shape) but would also be defined by the precise 
way stretch-injury takes place, i.e., the rate and extent of 
loading. Future studies may address this question in de- 
tail. 

Conclusions 

Neurites respond to sudden osmotic pressure changes 
with a fast volume regulation response. The initial phase 
is characterised by a nonlinear dependence of swelling 
rate on the initial osmotic pressure. Cytoskeletal per- 
turbation, especially microtubule disruption, accelerates 
swelling and increases the maximum volume reached, but 
does not affect the relaxation phase. Taking our results to- 



gether, we propose that mechanical forces due to the non- 
linear viscosity of the cytoskeleton slow down the initial 
phase of change of volume. This may provide instanta- 
neous integrity to the neurite while osmotic mechanisms 
"warm up". 
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